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Abstract 

Direct nuclear reactions d + A and p + (An) are described in the framework of 
three-body Faddeev-type equations. Differential cross section and analyzing 
powers are calculated using several optical potential models and compared 
with the experimental data. Quite satisfactory agreement is found except for 
few systematic discrepancies. 
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1. Introduction 

Direct nuclear reactions, dominated by three-body degrees of freedom, 
provide an important test for the models of nuclear dynamics. In the past, 
the reactions like the deuteron [d) scattering from a stable nucleus {A) 
were described using the approximate continuum discretized coupled chan- 
nels (CDCC) method Q. With few exceptions the spin degrees of freedom 
have been usually neglected. In such a case one can only calculate the un- 
polarized cross sections which are mostly sensitive to the central part of the 
nucleon-nucleus optical potential. However, there exist experimental data 
for polarization observables that depend as well on the spin-orbit part of the 
optical potential. The aim of the present work is to study the spin observ- 
ables m d + A and p + [An) reactions in the framework of exact three-body 
Faddeev/Alt, Grassberger, and Sandhas (AGS) equations 0, [sf such that 
all discrepancies with the experimental data can be attributed solely to the 
shortcomings of the used optical potential models. A part of the reactions 
to be considered here have already been studied by the CDCC method 
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Thus, the comparison with the results of Ref. [J] may allow us to draw some 
conclusions on the reliability of CDCC for calculating spin observables, in 
addition to the cross section benchmark [sj. 

The theoretical framework is shortly recalled in Sec. [2l the results are 
presented in Sec. [3l and the summary is given in Sec. HI 

2. Faddeev/AGS equations and dynamic input 

We describe d + A and p+ {An) type reactions using three-body {p, n, A) 
model. An exact treatment of the quantum three-body scattering problem is 
provided by both Faddeev 0] and AGS equations that are equivalent to 
the Schrodinger equation but are more suitable for the numerical solution due 
to the connectedness of the kernel. The Faddeev equations are formulated 
for the components of the wave function while the AGS equations, 



are a system of coupled integral equations for the transition operators Upa 
whose on-shell matrix elements (ipislUisalipa) are scattering amplitudes and 
therefore lead directly to the observables. In Eq. ([1]) Spa = 1 — Sf3a, Gq = 
{E + iO — Hq)~^ is the free resolvent, and T^- = Vfj + v^jGoT^ is the two-particle 
transition matrix, E being the available three-particle energy in the center 
of mass (cm.) system, Hq the free Hamiltonian, and Va the potential for 
the pair a in odd-man-out notation. The channel states for a = 1, 2, 3 
are the eigenstates of the corresponding channel Hamiltonian Her = Hq + 
with the energy eigenvalue E; thus, 1^^^) is a product of the bound state 
wave function for pair a and a plane wave with fixed on-shell momentum 
corresponding to the relative motion of particle a and pair a in the initial or 
final state. Observables of elastic scattering are calculated from the matrix 
elements with (3 = a while (3 ^ a corresponds to transfer reactions. 

We solve the AGS equations using momentum-space partial-wave basis 
where they become a system of integral equations with two continuous vari- 
ables, the values of Jacobi momenta. The employed numerical techniques are 
described in great detail in Refs. @, 0, @] • 

The AGS equations are applicable only to short-range potentials Va. Nev- 
ertheless, the long-range Coulomb force between charged particles can be in- 
cluded in this framework using the method of screening and renormalization 



3 




(1) 



cr=l 



2 



(gl-fiol. 11] which enables to calculate the Coulomb-distorted short-range part 
of the transition amplitude by solving the AGS equations with nuclear plus 
screened Coulomb potential; the convergence of the results with the screening 
radius has to be established. The method has been successfully applied to 



proton-deuteron [8|, |12| and a-deuteron [13[ elastic scattering and breakup, 
and to three-body nuclear reactions involving deuterons or one-neutron halo 
nuclei jsf. 

The dynamic input to the AGS equations are the potentials for the 
three pairs of particles. As the np interaction we take realistic CD Bonn 



potential [ij], in contrast to the usual CDCC calculations where a simple 
Gaussian np potential is used. For the nucleon-nucleus {NA) interaction, in 
order to study the model dependence, we use several different optical poten- 



tials, namely, those by Watson et al. |15[|, Menet et al. [16|, Becchetti and 



Greenlees [17[, and Koning and Delaroche [18|; the corresponding predictions 
in the following will be abbreviated by W, M, BG, and KD, respectively. 
Each of these potentials is fitted to the NA data in a limited mass and en- 
ergy range; those limitations are respected in the present calculations. The 
energy-dependent parameters of the potentials are taken at a half deuteron 
lab energy in the d + A reactions and at the proton lab energy for pA in 
the p + {An) reactions. In the latter case the nA potential is real and sup- 
ports a number of bound states corresponding to the ground and excited 
single-particle states of the [An] nucleus while all Pauli forbidden states are 
removed; the potential parameters and the resulting binding energies are 
given in Ref. [iQ] for ^'^C and ^''O nuclei. The interaction within np, nA, and 
pA pairs is included in partial waves with pair orbital angular momentum 
L < 3, 10, and 20, respectively, and the total angular momentum is J < 45; 
depending on the reaction some of these quantum numbers cutoffs can be 
safely chosen significantly lower, leading, nevertheless, to well converged re- 
sults. The pA channel is more demanding than the nA channel due to the 
screened Coulomb force, where the screening radius i? ?s 8 to 10 fm for the 
short-range part of the scattering amplitude is sufficient for the convergence. 

3. Results 

The experimental data for the spin observables are scarcer than for the 
spin-averaged cross sections. Nevertheless, a complete set of deuteron ana- 
lyzing powers over a wide mass and angular range is presented in Ref. 2^ 



for Ed = 56 MeV; other measurements exist as well but cover only rather 
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Figure 1: Differential cross section divided by Rutlierford cross section and deuteron 
analyzing powers for deuteron elastic scattering from ^^C and ^^O nuclei at Ed = 56 MeV 
as functions of the cm. scattering angle. Results for NA potentials W and M are given 
by dashed and dotted curves, respectively. The experimental data are from Ref. \2d\. 



narrow angular or mass range. We therefore concentrate first on 56 MeV 
polarized deuteron scattering from ^^C, ^^O, ^'^Si, ^°Ca, and ^^Ni nuclei. The 
results are presented in Figs. [1] and [21 




Figure 2: Differential cross section divided by Rutherford cross section and deuteron 
analyzing powers for deuteron elastic scattering from ^®Si, "'"Ca, and ^^Ni nuclei at Ed = 56 
MeV. Results obtained with NA potentials KD, M, and BG are given by dashed-dotted, 
dotted, and solid curves, respectively. The experimental data are from Ref. [21] for ^*Si 
and from Ref. [2^ for other nuclei. 
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As shown in Fig. [H in the case of ^^C and ^^O nuclei the predictions 
using the potential M describe the experimental data better than those of 
W, especially at larger scattering angles. There are significant differences also 
around Qc.m. = 22 deg where W predicts much deeper minimum in the cross 
section and thereby too sharp peaks in the analyzing powers; in contrast, the 
corresponding peaks obtained with M for ^^C are not sharp enough. Both 
potentials fail to reproduce the vector analyzing power Ay data at small 
angles Qc.m. < 40 deg; the reason for this failure, at least to some extent, may 
be the inability to describe Ay at small angles in p+A elastic scattering as can 
be seen in Ref . 15|] . The Ay and Ayy data at intermediate angles are slightly 



underpredicted, and, in addition. Ay for C shows a clear minimum around 
©cm. = 120 deg that is not present in the theoretical results. Furthermore, 
calculated A^z for both ^^C and ^^O nuclei at larger angles deviates quite 
significantly from the data. Otherwise the qualitative description of the data 
by the potential M is quite satisfactory. The comparison with the d+^'^C data 



from Ref. [22| at = 35 to 70 MeV that are limited to narrow angular range 
from 35 to 80 deg (not shown here) brings essentially the same conclusions. 

However, the potential M fails in the case of heavier nuclei ^'^Si, ^°Ca, 
and ^^Ni as shown in Fig. [2l It clearly overpredicts the differential cross 
section at Qc.m. > 40 deg while the analyzing power data are reproduced with 
similar quality as by other potentials. Quite surprisingly, at larger angles 
©cm. > 60 deg the old global potential BG describes deuteron scattering 
from ^°Ca and ^^Ni nuclei better than the new and precise KD potential 
which has many more parameters fitted individually to the nA and pA data 
of the considered nuclei. At smaller angles all potentials account for the data 
with comparable quality except for Ay where KD shows an additional sharp 
maximum instead of a minimum. The most serious discrepancies between 
theory and data take place in small- angle Ay and in large- angle Ar^z- While 
the latter one is similar as in the case of ^^C and ^^O nuclei, the former 
one is somehow different: the data are reproduced at very small angles, but 
then there is a narrow angular interval around Qc.m. = 30 deg where Ay is 
strongly overpredicted. Furthermore, in contrast to ^^C and ^^O, Ay in proton 
scattering from heavier nuclei is well described by the employed potentials. 
Thus, in this case the Ay discrepancy even partially cannot be explained by 
a poor NA data fit and must be due to another, presently unknown, reason. 
In addition. Ay and Ayy at larger angles are slightly underpredicted but the 
oscillating behaviour of the observables is reproduced rather well. Otherwise 
the qualitative description of the data, especially by the potential BG, is 
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Figure 3: Differential cross section divided by Rutherford cross section and deuteron vector 
analyzing power iTn for deuteron elastic scattering from ^^C and ^®Si nuclei at Ed = 29.5 
MeV calculated using the potentials W and KD, respectively. The experimental data are 
from Ref. 

quite satisfactory. 

Even larger discrepancy in the vector analyzing power Ay or, equivalently, 
iTii = Ay\^/2, is present at lower energies as shown in Fig. [3] for the elastic 
deuteron scattering from ^^C and ^^Si nuclei at = 29.5 MeV. Compared 
to Figs. [Hand [2], the discrepancy takes place at larger scattering angles, 30 to 
55 deg and 45 to 70 deg for ^^C and ^^Si, respectively. Outside those regions 
the calculations account for the data quite satisfactorily, at least in the case 
of ^^Si. The agreement for the differential cross section is rather similar to 
the one seen at Ed = 56 MeV. 

Deuteron scatteringfrom ^^O, ^°Ca, and ^^Ni nuclei at Ed = 56 MeV 
was calculated in Ref. [4|] in the CDCC framework using M and BG optical 
potentials. However, the np interaction in Ref. and in the present work 
is not exactly the same. In Ref. 0] the deuteron bound state was calculated 
with a realistic np potential but the np continuum was described by simple 
Gaussian potential acting in the triplet S and D waves only. Such a mixed 
potential cannot be included in consistent Faddeev-type calculations such 
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as of the present work. Nevertheless, since the considered observables are 
rather insensitive to the np potential, a comparison of the present results 
and those of Ref. jl] is meaningful. The agreement is good, there are only 
rather small differences, e.g., our Ay values for Qc.m. > 60 deg are higher by 
about 0.05 - 0.10 coming closer to the data. Thus, in addition to the cross 
section benchmark one can conclude that CDCC is quite reliable also in 
calculating the spin observables in deuteron-nucleus elastic scattering. 

The p+ (An) elastic scattering is a less interesting case since the differen- 
tial cross section and the proton analyzing power are quite strongly correlated 
with the corresponding observables in p + A elastic scattering to which the 
pA potential is fitted. The data for other observables like (An) analyzing 
power and spin correlation coefficients does not exist. Furthermore, the only 
available polarization data refer to proton-^^C scattering. Examples at pro- 
ton lab energy Ep = 17.5 and 35 MeV are shown in Fig. H] together with the 
differential cross section. The calculations describe the experimental data 
with the same quality as for the corresponding proton-^^C observables and 
even show similar deviations 1^, e.g., underestimation of the large angle 



cross section at 17.5 MeV and overestimation of the small angle Ay. How- 
ever, in contrast to the deuteron-^^C scattering, it seems that the potential 
W describes the data slightly better than M, except for forward angles. 

In the Faddeev/AGS framework the elastic scattering and the transfer 
reaction p + (An) ^ d + A are calculated simultaneously. Unfortunately, the 
data only exists for p+^^C — >■ d+^^C transfer at Ep = 65 MeV. Although this 
is slightly above the upper vahdity limit of the employed optical potentials 
W and M and both of them do not reproduce accurately the pA data at this 
energy, we show the results in Fig. [51 The description of the data is poor 
for both differential cross section and proton analyzing power, and is slightly 
worse when using the potential M. Note that a similar failure in the transfer 
cross section is observed in the d + p+^'^0 reaction at Ed = 63.2 MeV 



19 . 



4. Summary 

We performed Faddeev-type calculations of three-body direct nuclear re- 
actions in {p, n, A) model. The framework is the integral momentum-space 
AGS equations; the Coulomb interaction between the charged particles is 
included using the method of screening and renormalization; well converged 
results are obtained. Thus, all discrepancies with the experimental data can 
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Figure 4: Differential cross section divided by Rutherford cross section and proton analyz- 
ing power for proton-^'^C elastic scatte ring at Ep = 17.5 and 35 MeV. Curves as in Fig.[TJ 
The experimental data are from Refs. [2J| and [25| at 17.5 and 35 MeV, respectively. 
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Figure 5: Differential cross section and proton analyzing power for p + -'^^C ^ d + ^"^C 
reaction at Ep — 65 MeV. Curves as in Fig.[TJ The experimental data are from Ref. (26j . 
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be attributed to the shortcomings of the used optical potentials or even to 
the inadequacy of the three-body model. A realistic np potential and several 
parametrizations of the nucleon-nucleus optical potentials were used, all in- 
cluding spin-orbit interaction. Differential cross section and analyzing powers 
were calculated for the deuteron elastic scattering from ^^C, ^^O, ^^Si, ^°Ca, 
and ^^Ni nuclei, and for p + ^^C elastic scattering and transfer to d + ^"^C 
The description of the experimental data is mostly quite satisfactory, at least 
by some of the optical potential models, and is of similar quality as in the 
nucleon-nucleus elastic scattering to which the optical potentials are fitted. 
Systematic discrepancies are found in elastic proton and deuteron vector an- 
alyzing power Ay at small angles, and in the deuteron tensor analyzing power 
A^z at large angles. The former one, in the case of light nuclei ^^C and ^^O, 
may be related to a similar problem in the nucleon-nucleus Ay. The calcula- 
tions also fail in accounting for the available p + ^^C — > d+^^C transfer data. 
A satisfactory quahtative agreement with previous CDCC results ior d + A 
elastic observables is found indicating the reliability of the CDCC method 
for this type of reactions. 
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